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Electromagnetic field produced in high-energy heavy-ion collisions contains much useful infor- 
mation, because the field can be directly related to the motion of the matter in the whole stage 
of the reaction. One can divide the total electromagnetic field into three parts, i.e., the contribu- 
tions from the incident nuclei, non-participating nucleons and charged fluid, the latter consisting 
of strongly interacting hadrons or quarks. Parametrizing the space-time evolution of the charged 
fluid based on hydrodynamic model, we study the development of the electromagnetic field which 
accompanies the high-energy heavy-ion collisions. We found that the incident nuclei bring a rather 
^\ • strong electromagnetic field to the interaction region of hadrons or quarks over a few fm after the 

On I collision. On the other hand, the observed charged hadrons' spectra are mostly affected (Coulomb 

On ■ effect) by the field of the charged fluid. We compare the result of our model with experimental data 

and found that the model reproduces them well. The pion yield ratio -k~ /n'^ at a RHIC energy, 
$— ( , Au-|-Au 100+100 GeV/nucleon, is also predicted. 
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I. INTRODUCTION 
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C , Most of the important informations on the new phase of quark-gluon matter are obtained by investigating their 

v_H ■ space-time development in high-energy heavy-ion colhsions. These studies promise not only to specify the thermody- 
^^ [ namical conditions for the opening new phase but also to get useful knowledge of the evolution of the system. The 
^^ I Bose-Einstein correlations Q attract a great deal of attention as a powerful tool to study the geometry of the particle 
Qs [ emitting sources. However, very complicated interactions characteristic of these collisions prevent us from drawing a 
Qn^ . clear image of the sources. 

"^ ' The Coulomb interaction between a frozen out hadron and the strongly interacting charged hadronic matter may 

"^ , also be helpful in studying the evolution of the hadronic matter 0-0] . One of the advantages of studying its effects 

I ' is that there is no ambiguity in constructing the electromagnetic fields if one knows completely the motion of the 

O , hadronic or quark matter. Conversely, if one understands the electromagnetic fields precisely, one can know the 

2 ' behavior of the hadronic or quark matter in the high-energy heavy-ion collisions. Another advantage is that the 

electromagnetic fields carry information over the whole duration of the nuclear collision (actually, even of the stage 

before the collision.) Fortunately, a great number of charged hadrons are produced in the high-energy heavy-ion 

collisions. It is expected that one can obtain information on the electromagnetic fields by carefully studying their 

momentum distributions. The importance of the Coulomb effect at Bevalac energies has been pointed out by several 

C^ ^ authors ||-§. 

The hydrodynamic model |pO,nl| describes the space-time development of the hadronic matter and it can explain 
various experimental results. The purpose of this paper is to explore this model a little further into the problem of 
the electromagnetic fields due to the expanding charge distribution produced in high-energy heavy-ion collisions. 

In the following, we will confine ourselves only to symmetrical collisions such as Pb-|-Pb. The extension to non- 
symmetrical cases is straightforward. We will begin by describing in the next section how to compute the number of 
participant nucleons and their distribution in the transverse plane. In order to study the space-time development of 
the produced electromagnetic field, we shall assume a longitudinal scaling solution without transverse expansion for 
the charged fluid. We find this will be enough to reproduce the main features of the electromagnetic field in question 
and moreover, in computing the ratios of the yields (for instance 7r~/7r+), possible deviations are largely canceled out. 
In Section IE, we give account of how to estimate the net charge distribution inside the fluid along the collision axis. 
This will be done by using various experimental results on Pb+Pb 158 GeV/nucleon in the CERN-SPS. Then, in 
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Section IV, the electromagnetic field produced by this expanding fluid is computed. In addition to this, fields by both 
the incident nuclei and by the penetrating systems consisting of spectator nucleons are also taken into accountfj. The 
determination of the freeze-out hypersurface and computation of the transverse-momentum distribution is given in 
Section V. The results obtained are compared with the experimental data for the pionic yield ratio 7r~/7r+ by NA44 
Collaboration |12[ . Finally, we close our discussions with some concluding remarks and prediction for the pion yield 
ratio for Au+Au collisions at 100+100 GcV/nucleon. 

II. GLAUBER MODEL 

A. Participant and penetrating nucleons in A-B collisions 

The participant nucleon number in a nucleus- nucleus collision can be estimated by using the Glauber model [g_3|. 
To do this, let's first consider a h-A collision. At high energy, one can treat this problem as a collision of h with a 
tube in the nucleus (see Fig.l). Here, we assume that the cross section of the tube is given by the nucleon-nucleon 
inelastic cross section ain ~ S2mb. The probability that an inelastic interaction occurs in a volume clement <7inSz is 

SP{b, Zn) = CFinSz p{b, z„) = ain—p{b, z„) w -^ dz p{b, z) (1) 

where p(t>, z) is the nucleon density in the nucleus and TV is the number of equally divided volume elements in the 
tube (see Fig.l). The probability that h interacts with m elements in the tube is 

P{m; b) = ^C,„[JF(6)]"[1 - 5P{b)f-"' . 

Hence, we obtain for the total probability of occurrence of inelastic reactions in a h-A collision 

Pnel = Y. NCm[5P{b)r[l - (5P(b)]^-™ =!-[!- ^T{b) f (2) 

m— 1 

where T{b) is the so-called thickness function which is defined as T{b) = dz p{b,z). In the continuum limit, 
N —^ oo , this probability can be written as 

P,nei{b) = l-cxp[-a,„T{b)] . (3) 

As a natural extension of the ft,-tube collision, one can regard the nucleus-nucleus collision as a superposition of 
tube-tube collisions. Therefore, the total number of participant nucleons and penetrating nucleons originated from 
each incident nucleus is given respectively by 

Npartib) = Jdb' Tib') { 1 - exp [-araT{b' + b)]} , (4) 

Np,^{b) = Jdb' T{b')e^p[~a,^Tib' + b)] . (5) 

Assuming the nucleon density in the Woods-Saxon form 

Po 



p{b,z) 



1 + exp 



{^/z^ + b^-R)/d 



where d = 0.54 fm, Ra — tq x A^'^ and tq = 1.2 fm, the penetrating nucleon number in each hemisphere in Pb-|-Pb 
collisions turns out to be 30.7 and 58.4 for the impact parameter b < 3.5 fm (centrality 5%) and b < 6.1 fm (centrality 
15%), respectively. This means that the total number of participant nucleons in each case is N'P"-'^'^ =353.1 and 
j\fpart =297.7, respectively. See also Table H. It is gratifying to see that, when compared with the data (see Ng_g in 
the next section), these numbers are almost exact. 



^In a previous work by Barz ef aZ. H], it is not clear the role of the incident nuclei and penetrating systems. 



B. Transverse profile of the nucleon number density 

In the present work, we arc mostly interested in central collisions. So, we shall assume a cylindrically symmetrical 
distribution both of charge and of fluid density. However, in comparing the results with the data, we have also to 
take the experimental conditions into account. It is clear that the impact parameter is not always zero, but fluctuates 
from event to event both in modulus and direction. In order to simplify the computation, we shall symmetrize the 
fluid and charge densities by considering this fluctuation by averaging over b. Furthermore, we assume factorization 
of these densities with respect to the longitudinal and transverse variables. 

Let us estimate the transverse profile of this symmetrized nucleon number density. For a given nucleus- nucleus 
collision, whose geometry is shown in Fig. 2 (with a fixed 6), and, provided that the fluid contains nucleons which 
collide inelastically at least once, its nucleon density projected into the transverse plane (p = {x,y)) is obtained by 
adding 



7^(PT^)X 



l-cxp[-a,„r(p±-)] 



(6) 



So, by averaging over the impact parameter b of the incident nuclei 
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We also have, for the penetrating system and incident nuclei. 
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Observe that these quantities are related by 

Tfluidip) + Tpenip) = Tnucl{p) ■ 



(8) 



(9) 



(10) 



being the integral of the right-hand side equal to 2A. They are shown in Fig. 3. We can also define the normalized 
densities 



Hp) 



Up) 



dp p I dip Ti{p) 



(11) 



for the fluid, penetrating system and incident nuclei. In the following sections, we assume that the charge distribution 
in each system is proportional to these quantities. Also, the fluid density will be taken proportional to Tfiuidip)- 



III. NET-CHARGE DISTRIBUTION IN THE FLUID 



A. Space-time evolution of the fluid 



As mentioned in the Introduction, we approximate the evolution of the charged matter by a longitudinally expanding 
scaling solution, neglecting the transverse expansion. The transverse structure has already been discussed in the 
previous section. In such a case, the rapidity of the fluid is given by 



1, t + z 

■n — -in 

' 2 t-z 



It is convenient to take this and the proper time 



Vt^ 



together with p as the coordinates for describing the motion of our charged matter. Thus, 

t — T cosh 77, z = r sinh 77 , 



(12) 



(13) 



(14) 



for the time and the longitudinal coordinate in the ccnter-of-mass frame. Strictly speaking, the scaling solution is 
only valid if one has an "infinitely long" fluid (— cxd < t] < +cxd) with an infinitely large size and constant density in 
the transverse directions. For a real finite fluid, transverse expansion necessarily appears and also the longitudinal 
form suffers degradation beginning from the edges. In the present paper, we will neglect these effects, for the sake of 
simplicity because we find that they are not essential as far as the central region of the collision is concerned. 



B. Longitudinal profile of the charge distribution 



Let us now discuss the longitudinal structure of the charge distribution. The net charge distribution in the 77- variable 
(integrated charge distribution over p-space) in the fluid can be estimated in terms of the rapidity distributions of 
the produced particles 



dN„. 
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dN^- dN 
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dy 



dy 



(15) 



It should be noted here that the net charge distribution which we need to know is the one inside the fluid, so we have 
to estimate dNnet/dy soon after the freeze-out, i.e., r=10^ ^ 10^ fm. In this stage (before the hyperon decays), the 
multiplicity of the negative mesons {h~ = 7t~ ,K~) is slightly larger than that of the positive mesons (/i+ = 7r+, iir+) 
due to the large neutron-proton asymmetry in large stable nuclei whereas no such Coulomb-energy effect is expected 
in the expanding system. 

Consider the charge and the baryon number conservations in the collision: 
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where N^"-^'^ and TV^'"'' are the numbers of participant protons and neutrons, respectively, and y+,y° and Y~ stand 
for positive, neutral and negative hyperons. For the present, we shall concentrate on estimating the net charge 
distribution in the case of 158 GeV/nucleon Pb-|-Pb collision for the centrality 5%. From these two equations, we 
obtain the following formula for the negative charge excess in mesons: 
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(16) 



Here, we used the experimental results, Nn-n/Np-p~1.07, N^^_jt±/N/^_x ~ 0-3 Qill^l- In the last term above, 

T,^. One can then estimate the negative charge excess using the experimental value for the 

148.9 (See Fig.4) and Nj^_j^-- 



N^^-j^ contains also E° 



net baryon number Ng_g =352±12, being the net proton number Np^p = 



27.4 |15|, as 



N. 



-/i+= 



dy 



dNh- dN,,+ 



dy 



dy 
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(17) 



direct 



One can verify the consistency of this estimate, by computing the participant proton number iV|"""* which is approx- 
imately Z/A X Ng_g « 138.7 and comparing with the total net charge, obtained by using iVp_p data and eq. (Q^ 



iV, 



p-p 



N,-_ 



h+ 



138.9. 



Hence, by putting eq. ( p^ ) into eq. (|I5|) and by considering that the existing negative particle data include also p , 
finally we arrive at the following expression for the charge distribution of particles at the freeze-out time 
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(18) 



For Pb+Pb at 158 GeV/nucleon, we have N/i-j^p — 695±30 |1^ (see Fig. 5). The Np is estimated by using an empirical 
relation p/A= 0.12^^0.13 found in S* + S" and S + Ag collisions at 200GeV/nucleon ||ll,|l^. Because the multiplicity 
of A «7.5 (A/A «0.2) ||l|l, we obtain Na « 35 and Np « 4.4 (See Fig. 6). The solid line in Fig.7 shows the net charge 
distribution obtained by eq.(Qq), for the centrality 5%, in 158 GeV/nucleon Pb+Pb. 

Now, the rapidity distribution of produced particles given by eq.(lq) is a convolution of the charge rapidity distri- 
bution in the fluid and the thermal rapidity distribution of the produced particles 



dNnet .r f, dNnet 



1 



5{mT> cosh{y~ri)/T _ ^ ' 



(19) 



where T = 140 MeV is the freeze-out temperature and {ttit^) is the pion mean transverse mass, evaluated with the 
Bose distribution at the temperature T. So, in computing the electromagnetic field in the following section, we 
took a piece- wise flat rapidity distribution for the fluid net charge as seen in Fig.7 which, upon convolution e q.(|l9[ ), 
approximates as closely as possible the particle distribution of eq.(18). The normalization constant Af in eq.(|19|) is 
fixed by equating the integrals of dNnet/dy and dNnet/ drj. As for the fluid density, as mentioned in the beginning of 
this section we took a completely flat r/— distribution which gives the correct number of pions. 



C. Centrality dependence of the net charge distribution 

As is easily understood intuitively, the net charge distribution strongly depends on the centrality of the collision. In 
the previous subsection, we have estimated the net charge distribution by using the NA49 data on 158 GeV/nucleon 
Pb+Pb collision with the centrality 5 % jlq]. We would like to study the Coulomb effect on the tt^/tt^ ratio, which 
has been measured by the NA44 Collaboration |lj], considering events with 15 % centrality. The question is that there 
are no corresponding data which would allow us to estimate the net charge distribution under the same experimental 
conditions. We have mentioned in Section II and verified in Section 11 that the Glauber model we described there 
gives a very good account of the NA49 data. In view of this, we will account for the centrality dependence of the net 
charge distribution, by the use of the referred model. Concretely, by assuming that the longitudinal structure of the 
net charge distribution remains the same, changing only the transverse profile and consequently the normalization, 
we introduce the following factor to relating quantities with different centrality (x): 



dNj, 
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G(x') dNh 



centrality x' 



Gix) drj 



centrality x 



for h — h + p and p — p , where G(x) = N^""^*" for the centrality x ■ The participant nucleon-number dependence of 
hyperon multiplicity is studied in Ref. p9||. Using their results, we have 
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1.30 ,,r 
dNy 

drj 



centrality x 



for Y—A hyperon. The multiplicity values estimated for the centrality 15%, by using the data obtained for the 
centrality 5 % , are summarized in Table H. 

IV. ELECTROMAGNETIC FIELDS IN HIGH-ENERGY HEAVY-ION COLLISIONS 



Having been defined the charge distributions and their time development in the previous sections, 
compute the electromagnetic field produced by each one of such distributions. 



now we can 



A. Potentials of the longitudinal expanding charged matter 

An infinitesimal charge dQ moving along a trajectory r = ro(i) in the hadronic or quark matter (fluid) produces a 
Lienard-Wiechert four-vector potential 



A^'{t,x) 



dQ 



vdQ 



(20) 



,\x - ro{tO\ - v{t') ■ {x - ro{t')} ' \x - ro{t')\ - v{t') ■ {x - ro{t')}J ' 
where t' is the retarded time satisfying the condition, 

t' + \x - ro{t')\ =t . (21) 

In terms of the coordinates {T,p,ri), for ro{t') = (p', v^t') = (p', r'sinhyy'), eqs.(|2(]|) and (|2l| ) may be rewritten 



A^(^r,p,^) = ( ^°^^^^'^^ ,0, 0, '"^^'^' "^^ 



r cosh(?7 — 77') — r' ' ' ' r cosh(7] — 77') — r' 
where 

r'2 - 2t [cosh(?7 - 77')] t' + r^ - |p - p'p = . 
By solving eq.(E3h with the condition t > t' , we have 



t' = T cosh(77 — 77') — Y T^ sinh (77 — 77') + |p — p'P 

with T > |p — p'l . 

As discussed in the previous sections, the charge distribution in the expanding fluid is written as 

dQiv',p') rf. I n dN„et 

dp'drj' -^Tfi^^^^p)^-^, 

so the total potentials due to the whole fluid, A^ = 4>fiuid{T, p, ??) and A^ = Afiuid{T, p, 77), are given by 

, dNnet cosh 77' X 9{t -\p- p'l) 



(f>fiuid{T,p,'q)^ e dp' Tfiuidip') di] 



and 



respectively. 



Afiuid{T,p,ri)=e I dp' Tfiuidip') drj 



' a/t^ sinh (77 — 77') + |p — p'|2 
, dNnet sinh77' x e{T - |p - p'\) 



drj' 



T-^ sinh (7/ — 77') + |p — p'p 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



B. Potentials of the incident nuclei and the penetrating nuclear systems 

The incident nuclei (nucl) and the penetrating nuclear systems which are not included in the fluid (pen) also produce 
electromagnetic fields. Here, we assume that the penetrating systems consist of nucleons which suffered no inelastic 
collision and take their rapidities y the same as the beam rapidities ±yhe- The scalar and vector potentials can be 
obtained in a similar way as in the case of fluid. The potentials due to the incident nuclei are 



^nucl 



:r,p,r])^eZ fdp'fnuciip') Y. 

•^ n' = + , 



Jj' = ± Vba 



and 



Anuciir, p, 77)= eZ I dp' f^udip') ^ 



J)' = ± Vb^ 



cosh-q' X 9{\p — p'l — r) 
T^ sinh^(7; — 77') + |p — p'| 



sinhr/' x 9{\p — p'\ — r) 
T^ sinh'^(77 — 7;') + |p — p'l 



(28) 



(29) 



where Z is the proton number of each of the incident nuclei and T^uci^p') is the renormaUzed thickness function 
defined in eq.(|ll|). For the penetrating systems, we have 



and 



i>-pe-a{j, P, ??)= e {Z') dp' fpenip') ^ 

))' = ± Vbc 



ApeniT,p,Tj)^e{Z') dp' fpenip') Yl 
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cosh 77' X 9{t — Ip — p'\) 
r^ sinh {rj — rj') + \p — p'\'^ 



sinhr]' x 6{t — |p — p'|) 
r^ sinh (ry — r;') + \p — p'|^ 



(30) 



(31) 



where {Z') is the average number of protons which are included in each penetrating system and Tpen{p') is the 
renormalized thickness function defined in eq.dlT 



C. Total electromagnetic field in high-energy heavy-ion collisions 

The total scalar and vector potentials produced in high-energy heavy-ion collisions are now written as the following 
sums: 



<t>{'l',p,ri)= (f) fluid (t, P,f])+ <t>nucl (r, p, if} -f (j)pen {t, P, ??) , 

Az (r, p, 77) = Afiutd (t, p,ri)+ Anud {t, p,r]) + ^pen (t, P, tj) . 
Then, it is straightforward to obtain the corresponding electromagnetic fields by 

E{x,t) = -V(l){x,t) - —A{x,t) , 



where 



and they are explicitly given by 
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(37) 

(38) 

(39) 

(40) 
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(42) 



D. Space-time structure of the electromagnetic field 

In order to visualize the space-time distribution of the electromagnetic field in high-energy heavy- ion collisions, 
we show now some results of numerical calculations by means of eqs.(36)-(41). First, Fig. 8 shows the transverse 

component of the electric field, Ep (= \ E"^ + E^), as function of p (= y/x^ + y^) at t = 15 fm in the center of mass 

system. One can see a high peak in the electric field which increases a.s z ^ t, i. e., as the proper time r — > 0. This 
peak is produced by one of the incident nuclei before the collision. For smaller values of z, the peak which appears 
close to p ^ Rpb (radius of Pb) is produced by the expanding charged fluid. The time dependence of -Ep at p = 5 fm 
for various z are shown in Fig. 9. One can observe the propagation of a sharp peak produced by one of the incident 
nuclei. The apparent discrepancy of the result for z = is due to the fact that both nuclei contribute to the field in 
this case. It is also found that the width of the peak is about At ^ 3 fm. 

The longitudinal component of the electric field E^ is shown in Fig. 10. The magnitude of E^ is much smaller than 
that of Ep. At z—Q, E^ =0 because of the symmetry. Except for the case of z = 14 fm at small p, the sign of E^ is 
negative mainly due to the charge concentration at large- ry region (see Fig. 7). The positive Ez in z = 14 fm is due 
to one of the incident nuclei. 

The magnetic field at t = 15 fm is shown in Fig. 11. Because of the expansion of a positively charged matter, the 
sign of B^ at z > is positive. One can also roughly estimate the maximum value of the electromagnetic energy 
density Emax- The maximum value of Ep « 7.5 [MeV/fm] is reached on z ~ fm and p ^ Rpb at i ~ 1 fm. Then we 
have e„iax ~ 0.3 MeV/fm^. 

V. FREEZE-OUT SURFACE 

To obtain the freeze-out surface, we use the particle density in the fluid parametrized as described in Section IE, 

n(r,ry,p) = f/,„,(p)^5^xi, (43) 

di] T 

where dNfiuid/drj is constant in some rapidity interval. Since the most of the produced particles are pions, let us 
reinterpret it as the portion which corresponds to pions and fix the parameters by fitting the data on 

^«3x^, (44) 

ay dy 

where dN,r/dy is related with dNfiuid/dij through 

dN^ r f dNfiuid 1 ..(.x 

-^-jy J dri—^ X ^^^^^ ,,,h(,-,)/T _ 1 • i45j 

Here we are including only pions, considering the other hadrons as contaminations. The NA49's data for negative 
hadrons include p and K^ |13]. Using the experimental results in the mid rapidity region, {K^ + K~ -t- 2K^) / {-n^ + 
tt" + TT^) « 0.135, K^ /K^ Ril.8 and K^jK^ Ril.8 [Q, one can estimate the multiplicity of kaons as shown in Table. 
H and the total pion multiplicity (including vr") is 1950 and 1640 for the centrality 5% and 15%, respectively. The 
solid line in Fig. 12 shows the total pion rapidity distribution for centrality 15%. The shaded area and the dashed 
curve are our parametrization for dN ji^idj dr\ and dN^^/dy respectively. 

Now, we can compute the freeze out surface by equating the density given in eq.(|43|) to the critical value Uf for the 
Bose-Einstein distribution at the freeze-out temperature T/ 
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T ■' 2tt^ 
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For r/=140 MeV, one obtains rif— 0.096 fm ■^. The intercept of the freeze-out surface in (r, p)-plane is shown in 
Fig.l3. 



VI. SINGLE-PARTICLE SPECTRA AND COULOMB EFFECT 

A. Transverse-momentum spectra of charged particles 

Now that the electromagnetic field and the freeze-out surface have been obtained, one can solve the equation of 
motion for a classical charged particle, 



mv{t) 



dty/l-v^{t) 



q[Eix,t) +v{t)x B{x,t)] 



(47) 



starting from the freeze-out surface. In the present work, we have neglected the effect of the Coulomb field before the 
freeze out. The time evolution of the tt"^- and the 7r~- momenta, as given by eg. (|4^) , can generally be expressed in 
terms of the location and the initial momentum rriTO at the freeze out as 



mj = m^{p,ri;mT,o,yQ,^), 



(48) 



where p,ri,mrj.Q,yQ,ip are the phase space coordinates of tt at the freeze-out time. We show in Fig. 14 some results 
for the specified initial conditions. One can see that tt^s are accelerated outward and tt^s are decelerated by the 
electromagnetic field. The rapidity j/o also shifts to y± . However, this change is quite small, \y± — ya\ ^ 0.1 . 

Since the particle number density is given by eq.(B3), if we neglect the Coulomb effect, the invariant cross section 
for TT '*^ would be written in the form 



E 



1 



atot <Pp ~ (27r)3 
1 



(27r)3 



P^ f^'^ ^ emT0COsh{ri-y)/Tf _ I 

F.O. surface ^ ^ 

drj < pdp Tf{p) mTo cosh(77 — y) + I dTTpf{T) Pto cos{(p 

Vmin I J 



1 



So, the transverse momentum distribution in this case is 



dN 



niTodm-j 



(2^)2 



dydip / dr] < pdp Tf{p) ttIto cosh(ry — y) -f dTTpf{T) pT^oCOSip 



amTOCosh{ri~y)/Tf _ ^ 



omro cosh(?)— j/)/T/ _ ^ 



Since the transverse mass ttito at freeze-out time shifts to m^ as given by eq. ([l8[) , the single particle spectrum 
eq.dSQ) is distorted for 7r+ as 



diV^ 



mtdmi 



dydcp / dry pdp Tf{p) 



(27r) .. 

X drriT-o S[mT,o-m+{p,?],m-r,y,ip)]^ m-ro coah{r] - y) + p^o 



dTf{p) 



1 



{2ny 



dydf / drj pdp Tf{p) 



^m^ 



dmi 



{ m^ cosh(77 — y) + p.^ 



dp 



cosyj 



1 



mTO=™T 



dp 



cosip } 



^m-To cas\i(ri -y)/Tf _ ^ 
1 



(49) 



(50) 



gm+ cosh(r;-y)/T/ 



(51) 



For TT^, one can obtain a similar expression. 

However, these are not yet the final results for the charged-particle spectra. In writing eq. (pl|) , we have neglected 
the possibility of frozen-out particles be reabsorbed by the fluid again. So, we have checked all the trajectories of 
TT^ mesons after leaving the freeze-out surface and excluded those ones which crossed it in, considering them as 
re-absorbed by the fluid. 



B. Yield ratio tt /tt^ 

Using the final results of eq.(pl|) which we have just derived, we obtain a formula for the pionic yield ratio 



■n /tt^ 



1 dN- 



1 dN+ 



rrirj. dnirj. / mt dmif 



(52) 



To compare more precisely our results with the NA44 data |T^ for centrality 15% , we introduce the acceptance 
correction factor A{mrj,,y) to our formula eq.(|5l|), which would approximately reproduce the nature of the NA44's 
detectors: 

A{m.^;y) = 6{y + l.67^/ml -ml- 4.2) . (53) 

Here we neglect the rapidity shift due to the Coulomb field. 

The results of our calculation together with NA44's data are shown in Fig. 15. As seen, our results approximately 
reproduce the data. We shall now look more carefully into our results in Fig. 15. It is found that our model overestimate 
in the region m-Y-m-,: ^ 5MeV and underestimate in the region 10 < mT-rriT^ < 80 MeV. The underestimation of our 
model may be improved by introducing a transverse expansion for the charged fluid, because it makes the freeze-out 
time shorter and then the frozen out particles would feel a stronger electromagnetic field. Another factor which we 
have neglected in the present calculation, namely the effect of the field before the freeze out, also contributes to correct 
the results in the right direction. It has also been suggested by other authors p^ that 7r^/7r+ ratio has origin in 
hyperon decays f\ 

Save particle creation and annihilation, no quantum effect has been considered in the present calculation. Due to 
the uncertainty relation, evidently the classical calculation only has meaning if the wave packets are sufficiently small. 
This implies that the Coulomb effect for small-px particles [§|,^ would be smaller, because such particles would have 
a large wave packet for a good definition of the momentum and would feel a weaker field, which is an average field 
over their extension. This deviation would be larger for n^ than for 7r+, because the former is attracted toward the 
symmetry axis whereas the latter is repelled away. 

We also show our result of 7r^/7r+ for the centrality 5% in Fig. 16. In this calculation, the rapidity of the pions is 
fixed at y,r—0-0 and 2.0 and set A{mrj,; y) = 1. 

VII. CONCLUDING REMARKS AND DISCUSSIONS 

In this paper, we have investigated the structure of the electromagnetic fields which accompany high-energy heavy- 
ion collisions and computed the Coulomb effect for pion spectra caused as a consequence. We have considered three 
possible origins of the total electromagnetic field, i. e., i) incident nuclei, ii) longitudinally expanding charged matter 
and iii) penetrating systems which are composed of surviving protons. Glauber model has been used to estimate the 
number of participant nucleons and to determine the transverse profile of the charged system after the collision. The 
centrality dependence of the net-charge is also taken into account. In the description of the space-time evolution of the 
charged matter, a parametrization based on one-dimensional hydrodynamical model has been used, with a constant 
rapidity rj. 

The main results of this paper could be summarized as the following. 



• 



• 



Compared with the expanding charged matter, the incident nuclei bring a rather strong electromagnetic field 
into the interaction region in r '^ a few fm. 

The NA44 Pb-|-Pb data on 7r^/7r+ are approximately reproduced by our model (especially due to the longitudi- 
nally expanding fluid part). This means that the longitudinally expanding charged matter plays an important 
role on the final-state interaction of produced hadrons. 

• As expected, the penetrating systems seem to be not important for the observables in the mid-rapidity region. 
(However, it should be noted that they could have non-negligible contributions in the fragmentation regions.) 

In addition to these points, we would like to note that the Coulomb effect seen in the yield ratio 7r^/7r+ may 
depend both on the dumping speed of the electromagnetic field and on the cooling speed of the hadronic matter. 
This point leads us further into consideration of the Coulomb effect which could be observed in the RHIC energies. 
One would expect that the strength of Coulomb effect may be small, because the fluid takes much time to cool 
down before the freeze out begins and the net-charge density in mid rapidity region become small at that time. 



■^Although most of vr^ from A-decay should appear as small p-r particles, in the NA44 Collaboration's setup almost all the 
hyperons (with rapidity in the range yiab. < 2~3) decay in the first dipole magnet region |21|. This means that very few 
daughter tt~ can enter their detector. So, although we think that this effect should be studied, it is seemingly very small. 
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To see this quantitatively, we present the prediction for the Coulomb effect in Au+Au collisions (centrality 5%) at 
lOOGeV+lOOGeV energy in Fig. 19. The net-charge distribution and the fluid distribution used in this calculation 
are shown in Fig. 17 and Fig. 18, respectively. For the net-charge distribution, we assume that the total net charge 
and width of 'shoulders' at large jryj regions is to be the same as that of SPS case^.The total number of pions in the 
RHIC energy in Fig. 18 is determined by the energy conservation law for the total system. As shown in Fig. 19, the 
enhancement of 7r^/7r+ is seen at very small ttit -Ti^r region only. 

As remarked in Section VI, there remained some open questions: effect of the transverse expansion, charge redistri- 
bution inside the fluid, introduction of the resonance decays and the effect of the uncertainty relation. These questions 
are presently under investigation. 
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Pb+Pb 158 GeV/nucleon 



P-P 

h~ = n~ 

A-A 

A 

P 



5 % centrality (fe <3.5 fm) 



15% centrality {b <6.1 fm) 



.K- 



352 

148.9 

695 

27.4 

7.5 
4.4 



iVs-Rin Ref. ^ 

Ref. Ili 

Ref. [15] 

Ref. jij] 

Ref i|iJl 

empirical 
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125.5 

585.9 

22.5 

5.6 

3.5 



Glauber 

Glauber 

Glauber 

Ref. Ili + Glauber 

Ref. 119+ Glauber 

Ref |19 +Glauber+empirical 



N„et = 2Z 

<Z' > 
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12.5 



single hemisphere 



10.6 
116 
24.6 



single hemisphere 



K- 



81 
45 



68 
38 



Glauber 
Glauber 



TABLE I. 

Table caption 

Table I: Centrality dependence of the participant nucleon numbers, hadron multiplicities and net charge excess in 
Pb+Pb collision at 158 GcV/nucleon. 
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Figure captions 



Fig.l 
Fig.2 
Fig.3 



h-tuhc interaction in high-energy h-A colhsions. 

A-B colhsion with impact parameter b as seen in the transverse (x-y) plane. 



Transverse profile of the average nucleon density for the fluid, penetrating systems and incident nuclei, 
respectively. The results are shown for Pb+Pb collisions (ccntrality 5%). 



Fig.4 
Fig.5 
Fig.6 
Fig.7 



Net proton {p-p) distribution in Pb+Pb 158GeV/nucleon. 

Negative hadron (tt^, K^ and p) distribution in Pb+Pb 158GeV/nucleon. 

A distribution and estimated p distribution in Pb+Pb 158 GeV/nucleon. 



Net-charge distribution of the fluid (shaded area). The dashed curve is the result of convolution, eq.(19). 
The solid curve is the charged particle distribution, eq.(18). 

Fig.8: Electric fleld Ep at t=15fm. 

Fig. 9: Time dependence of Ep at p~5ini. 

Fig. 10: Longitudinal component of the electric fleld E^ at i=15fm. 

Fig.ll: Magnetic fleld B^ at i=15fm. 

Fig. 12: Total pion distribution and fluid distribution dN fluid/ drj. 

Fig. 13: Freeze-out surface projected on the T-p plane. Due to our parametrization, eq.(|l2[), the curve in r-p plane 
is independent of rj. 

Fig. 14: Time evolution of the momentum of 7r+ (solid line) and 7r~ (dashed line). Initial conditions are chosen as 
a;=6fm, y—z—Qiva at i = 15 fm with momentum, 20, 40, 60 and 80MeV in the ~\-x direction(transverse radial 
direction). 

Fig. 15: Comparison of our result for the yield ratio tx^ jix^ (solid line) with NA44 CoUab's data |12|. Remark that 
we did not flt the data, but the normalization of the curve has been done by estimating the overall ratio through 
data. 

Fig. 16: Our prediction of 7r^/7r+ in Pb+Pb 158GeV/nucleon for the centrality 5%. The pion rapidity y^r is flxed 
to be 0.0 and 2.0. 

Fig. 17: Net charge distribution at a RHIC energy 

Fig. 18: Fluid distribution at a RHIC energy. The pion multiplicity is determined by the energy conservation law 
for the total system. 

Fig. 19: Prediction for 7r^/7r+ at a RHIC energy. 
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Pb+Pb Collision (centrality 5%) 




p[fm] 



10.0 



FIG. 3. 



50 



Pb+Pb 1 58GeV/nucleon (centrality 5%) 



40 






NA49 Collaboration 
p-p= 148.9 




0.0 

rapidity y 



full phase 
\ =153.0 



4.0 



FIG. 4. 



15 



200 



T3 



100 



Pb+Pb 158GeV/nucleon (centrality 5%) 



NA49 Collaboration 




0.0 
rapidity y 



FIG. 5. 



T3 



8.0 


Pb+Pb 1 58GeV/nucleon (centrality 5% ) 


NA49 Collabolation 


6.0 


n r 






n 


^A=7.5 


4.0 


F 


- - 


\ 


2.0 


f'\ 


n n 


1 \ ^^ P=4.4 

J V 



-4.0 



-2.0 



0.0 
rapidity y 



2.0 



4.0 



FIG. 6. 



16 



T3 







Pb+Pb 158GeV/nucleon (centrality 15%) 




40 






' 1 ' 1 
eq.(18) 










eq.(19) 


30 


/ 


\ 


1.4 

ilillliili 


i jiiHiii 


20 


O 






\/ 


\ 




CT) 


1 




V '^ 




10 


- 






\ CD 


- 


n 


\ 


J 




V 


/ 



-4.0 



-2.0 



0.0 
11 



2.0 



4.0 



FIG. 7. 



Pb+Pb 1 58 GeV/nucleon (centrality 5%) 



z=14fm 




50 100 

P [fm] 



150 



FIG. 



17 



Pb+Pb 158GeV/nucleon (centrality 5%) 




0.4 



0.2 



FIG. 9. 



Pb+Pb 158GeV/nucleon (centrality 5%) 



> 



N 
LU 



0.0 



-0.2 



z=Ofm 



-0.4 




10 20 



30 40 50 60 70 
P [fm] 



FIG. 10. 



18 



Pb+Pb 158GeV/nucleon (centrality 5%) 



> 

CD 



DO 




50.0 100.0 

P[fm] 

FIG. 11. 



150.0 



600 
500 
400 



Pb+Pb 158GeV/nucleon (centrality 15%) 



— pions observed 
pions from the fluid 



fluid 




FIG. 12. 



19 



Pb+Pb 158 GeV/nucleon (centrality 15%) 




0.0 



2.0 4.0 6.0 

p[fm] 



8.0 10.0 



FIG. 13. 



100 



freeze out point: x=6fm, y=z=Ofm, t=15fm 



o 
> 



80 



60 



initial px= 80MeV 



initial px= eOlVleV 



Q. 



40 



20 



initial px= 40MeV 



initial px= 20MeV 



% 



rf 







100 200 

t [fm] 



300 



FIG. 14. 



20 



2.0 



1.8 



1.6 



+ 


1.4 


'c. 


1.2 




1.0 




0.8 



0.6 




200 300 400 

m^-m^ [MeV] 

FIG. 15. 



Pb+Pb 158GeV/nucleon (centrality 5%) 



3.5 



+ 2.5 



1.5 



y=0.0 

■'71 




y=2.0 ^ ^***C.-.o«. 



0.5 



10 20 30 

m^- mjj[MeV] 

FIG. 16. 



40 



50 



21 



Au+Au 100+100 GeV/nucleon (centrality 5%) 
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